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The nuclear incompressibility re is investigated in asymmetric systems in a mean field model. The 
calculations are done at zero and finite temperatures and include surface, Coulomb and symmetry 
energy terms for several equations of state. Also considered is the behavior of the incompressibility 
at constant entropy kq which is shown to have a very different behavior than the isothermal re. 
Namely, kq decreases with increasing entropy while the isothermal re increases with increasing T. 
A duality is found between the adiabatic kq and the T = isothermal re. Analytic and also simple 
approximate expressions for re are given. 



I. INTRODUCTION 



Heavy ion collisions produce hot and dense hadronic matter [l|, |2| . Such matter is also encountered in nuclear 
astrophysics and in supernovae explosions. Future rare isotope accelerators will explore nuclei at the limit of isospin 
asymmetry. In neutron stars one also encounters system with a large neutron to proton ratio. Thermodynamic 
quantities play an important role in characterizing properties of such systems. One such quantity is the nuclear 
incompressibility. While the nuclear incompressibility at zero temperature has been studied for an extended period 
1111 it is only relatively recently that its temperature dependence has been of concern 0, ■ In this paper 
we will focus on the temperature dependence of the incompressibility in asymmetric systems, including the role of 
surface, symmetry and Coulomb terms. We will also study the adiabatic incompressibility which at T = would be 
the same as the isothermal incompressibility, but differs from it at finite T. The isothermal incompressibility will 
increase with T while the adiabatic incompressibility will decrease with entropy |2j . 

As a baseline, we will begin with a mean field discussion of its behavior with T to see the contributions to incom- 
pressibility of effective mass, excluded volume, symmetry energy, surface energy and Coulomb energy. In Sect. II, we 
will summarize a general description of incompressiblity and results of Ref.0 which considered the temperature and 
entropy dependence of isothermal and adiabatic incompressibility of infinite symmetric nuclear matter with a simple 
Skyrme interaction. In Sect. Ill and IV we will consider the contributions of effective mass and excluded volume. 
In Sect.V we consider the incomprssibility of asymmetric finite nuclear system which has surface energy, symmetric 
energy and Coulomb energy contribution. The conclusion is given in Sect. VI. 



II. GENERAL CONSIDERATIONS 



First, we define a quantity k, incompressibility coefficient, as 



a d 2 (E/A) 2 d 2 (E/A) n , r2 d\E/A) u2 d 2 (E/A) 
This quantity is evaluated at the saturation density po where E/A has a minimum. The giant monopole resonance 



energy is related to re and this relation is Eq = y T^py • If the temperature is kept constant in the above derivatives 

we have the isothermal incompressibility re, and if the entropy is held fixed, the result is the adiabatic incompressibility 
kq. This re is related to the second order Taylor expansion coefficient of E/A with p or E/A = a + b(p — pq) 2 /2. 
Thus re = 9pyb. The quantity re defined above is not the isothermal compressibility defined in thermal physics as 

K = -(±)(dV/dP) T with T held fixed and here P is the pressure. Since P = -dF/dV\ T , we have K = (i) ■ 

At T = 0, F = E — TS = E, and thus re = 9/(poK). This reciprocal connection between K and re is no longer true 
at finite T. A reciprocal relation would exist between K and a re defined by an equation similar to Eq.lQ, but with 
energy E replaced with the free energy F. The K can be connected to the density fluctuations through the equation 
(p 2 ) — (p) 2 = ({p) 2 /V)TK Near a first order phase transition and at its critical point, the K can become very large 
since small changes in the pressure can produce large changes in the volume. Large density fluctuations at a critical 
point give rise to the phenomena of critical opalescence where the presence of droplets of all sizes scatter light. In the 
limit of an ideal gas K = 1/P. For a Van der Waals gas, whose equation of state is (P + a(N/V) 2 )(V — bN) = NT, the 
isothermal compressibility K will diverge above and near the critical temperature T c as 1/(T — T c ) 7 , where the critical 
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exponent 7 = 1 which is the mean field value. Its reciprocal will go to zero. The energy per particle for the Van der 
Waals gas is E/N = (3/2)T-a(N/V) for V > bN, while the entropy per particle is S/N = log[(V ~bN)/ (NX 3 )} + 5/2. 
The A is the De Broglie or quantum wavelength A = h/(2irmT) 1 ^ 2 . The energy per particle monotonically decreases 
with increasing density p = N/V and is at it's minimum value at the smallest density determined by the excluded 
volume V = bN. Besides the isothermal compressibility, an adiabatic compressibility can be obtained by keeping the 
entropy constant. The adiabatic compressibility is Kq = — (y)(dV/dP)s with the derivative at constant entropy. 
The reciprocal is related to kq, the adiabatic incompressibility as kq = 9/(pqKq). The energy per particle of a Van 
der Waals gas can have a minimum with V at constant entropy. Since the natural variables for energy are entropy 
and volume from dE — TdS — PdV variations of the energy with V at constant S bear a similar relation to variations 
of the Hclmholtz free energy with V at constant T where dF = —SdT — PdV . Thus the adiabatic incompressibility 
of Eq.m can go to zero for a Skyrme interaction as we shall see. The minimum point (also maximum point) in the 
energy occurs at zero pressure since dE/dV at constant entropy is —P. Therefore E at constant S has the same 
maximum and minimum points with variations in V or R or density as F at constant T since both derivatives are 
— P which is set to 0. 



A. Review of Infinite Matter Incompressibility 

Our mean field discussion is based on a Skyrme interaction. The Skyrme interaction shares some common features 
with the above Van der Waals equation such as an intermediate density attraction (the a(N/V) 2 term above) and a 
higher density repulsion (the excluded volume term in Van der Waals). The excluded volume leads to a much stronger 
repulsive term than the Skyrme repulsive term since it involves an expansion to all orders in the density. To keep the 
discussion simple to begin with, we consider uncharged symmetric nuclear matter with no surface energy terms. The 
Skyrme interaction energy is then 

U/A= -a Q p + a aP 1+a . (2) 

The ao term gives a medium range attraction while the a a term is a short range repulsion. 

At T = 0, the kinetic energy Ek/A — (^)Ep(p). The coefficients ao and a a are fixed to give a binding energy per 

particle E B /A = 16MeV at saturation density p = p = 0.15/fm 3 . At this density E F = ^;(^p) 2/3 = 35MeV and 
E K /A = 21MeV. From E = E K + U and d(E/A)/dp = at p Q the coefficients are 

a Po = b Q = (E K + E B )/A + b a = {E K + E B )/A + {l/a){E K /3 + E B )/A = 37 + 23/a 

aa a pl +a = ab a = (E K /3 + E B )/A = 23. (3) 
2 

(l + a)b a = b --E K /A 

The incompressibility coefficient k at T = is then 

K = -2E K /A + 9(1 + a)aa a pl +a = E K /A + 9E B /A + a(3E K /A + 9E B /A) = 165 + 207a. (4) 

For a — 1/3, k — 234MeV at po- Smaller values of a lead to softer equations of state and lower n. In the limit 
a — > 0, logarithmic terms appear in Eq.|J5J| coming from the presence of a factor (x/a)(l — x a ) — (x/a)(l — e alogx ) — > 
-xlog(a;). The x = p/p = (Ro/R) 3 with R% = A/(p 4n/3). The a -> limit is the softest EOS allowed by Eq.Q, 
and this limit gives from Eq.QJ a value of k = 165 MeV. A stiff EOS has a = 1 and k — 372 MeV. Recent calculations 
|9llTollTlj| have a value of k = 210 ~ 270 MeV and suggest a value of a = 1/3 in a Skyrme type approach. A larger 
range of values of k, from 211 ~ 350 MeV, were reported in Ref.^l- Because of these uncertainties in n from more 
realistic forces, we will present results for various values of a, from a w to a = 1. 

B. Corrections to the Infinite Matter Result 

If the corrections to the infinite nuclear matter incompressibility from finite size terms or non-zero temperature or 
entropy terms are small, then these corrections can be obtained by using the following method. Let Eq(R) be the 
nuclear matter energy per particle EOS which has a minimum at Rq and an incompressibility kq. If we add to this 
a term E X (R), so that E(R) = Eq(R) + E X (R) then the minimum shifts to a new point R m = Rq + AR X . The new 
minimum and K can be found by making a Taylor expansions around of Rq. The shift is given by 



AR x /Ro = -R (dE x {R)/dR)/K 



(5) 
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with the derivative evaluated at R since R 2 d ^/ 2 A \ft m ~ R 2 d dtt^ \iio f° r sma U Aiij.. The new k is 

S4 = ^(WD 



di? 3 



^ d 3 (E (R)/A) 



The various quantities are evaluated at Rq and each ' represents one derivative wrt R. Corrections to k involving 
skewness Sk or third derivative of the energy were pointed out in Ref. [4|, |5j . However note here that the skweness 
Sk is defined with opposite sign from others. Ellis et al used the correlation between compression modulus 
and skewness coefficient to examin the implications in a relativistic Hartree-Fock approximation where the E x is the 
Coulomb interaction. The above expression is a modified version of their result. When two small correction terms 
are present then the above result is easily generalized and each adds its own separate contribution. If the corrections 
are large then the above approximate expression will not be very accurate and it is best to calculate the shift and k 
exactly. We will use these results for the shift and k to find simple expressions for k for the general case of arbitrary 
a. Eq.Q) gives an expression for the incompressibility at T — in terms of Ek and Eb- This will be kq = Ko(a) in 
Eq.0 and Eq.lJBJ. The skewness Sk = R 3 (E (R)/A)"' , in terms of E K , E B and a, is 

Sk{a) = -27(E B /A)(3 + 4a + a 2 ) - {E K /A){ll + 36a + 9a 2 ) = -3(509 + 828a + 207a 2 ). (8) 

The last equality in Eq.JSJ) follows when Ek = 21MeV and Eb = 16MeV. Comparing this last equation with the k 
of Eq.QJ we see that the ratio of Sk/ k is of the order of 10 and somewhat insensitive to a. In particular for a = 0, 
1/3 and 1 this ratio is —9.25, —10.36 and —12.45. We will apply this procedure or exact solutions to evaluate the 
contribution of surface, Coulomb, symmetry energy, temperature and entropy to the nuclear incompressibility. For 
the particular case a = 1/3, exact solutions are also easy to obtain in several situations which lead to a quadratic 
equation for the equilibrium point R when E is written in terms of R. 



C. Non-zero Temperature and Entropy Considerations 

At non-zero T and with T < E F , the kinetic energy is E K /A = (^)E F + (^-)T 2 /E F H The energy per 
particle is 

El A = 21^/3 + |!_ _ bQX + baX i +a . (9) 

This equation can not be used for very small x since the condition T <C Ep will fail. We can use the procedure of 
Sect|ni20- Taking E X {R,T) = 0.0517T 2 i? 2 /A 2 / 3 in Eq.® and Q, we obtain the following expression for k and 



K(a) = R 



2 d 2 (E(R)/A) 



dR 2 



= K (a) - 0.1034T 2 i? 2 /A 2/3 {S*fc(a)/K (a) + 1} (10) 

T 

2/3N, 



x m {a) = Pm /p = (Ro/R m y = (1 + AR x /Ro)- d = (1 - O.IQMT'/A 2 ' 3 )- 3 (11) 

Note here that x m > 1, i.e., the saturation density is larger than p in an isothermal process. Since Aj (AttRq/'S) — 
po = 0.15/fm 3 and thus R 2 / 'A 2 ! 3 = 0.734fm 2 we have the following final simple results: 

n(a = 0) = 165+ 1.16T 2 
n(a = 1/3) = 234+ 1.32T 2 

K(a = 1) = 373+ 1.61T 2 (12) 

From Ea. (|12j) we see that the first term is very sensitive to a but the finite temperature correction is somewhat 
insensitive to a pj. Since Sk/n® is insensitive to a with a value of about —10, the second term does not scale with 
the first term in Ea.l|T2"fl as a changes. At T — 7.5MeV the T dependent corrections in Ea. lfT2^l are 65, 74 and 91MeV 
while the first terms are 165, 234 and 373 MeV for a = 0, a = 1/3 and a = 1. These results are very close to the 
exact solution result of Ref. p| • The exact result is given by 

k{T) = -42x 2 { 3 + 0.705T 2 /x 2 { 3 + 9a(l + a)b a x]+ a (13) 
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At T = 2.5, 5, and 7.5MeV, and for a = 1/3, the values of re are 242, 265 and 302MeV, respectively. The value of 
a = 1/3 and T — is also useful. It leads to an energy whose minimum can be obtained from a quadratic equation 
in the radius using p = A/(4nR 3 /3). Namely E(R)/A is 

/ R \ 2 f R \ 3 /f?\ 4 4 2 / 3 A 4 4 / 3 

E(R)/A = 21 (^j - 106 (^J + 69 (lfJ = 28.626-^-168.7-^ + 128.2-^ (14) 

At finite T a term (7r 2 /140)T 2 (i?/i?0) 2 = 0.0517T 2 i? 2 /A 2 / 3 is added to the right side of Eq.©. When T is replaced 
with entropy per particle 5/ A then this T dependent term becomes 4.836(S , /A) 2 A 2 / 3 /i? 2 since S = (tt 2 /2)TA/E f at 
low T. This S/A term can simply be added to the first term on the right side of Eg. (11411 since both have the same 
A 2 / 3 /R 2 dependence. In the limit a — > 0, x m satisfies the equation 14(a; 2/3 - x) + 23x log x = (it 2 /210)T 2 /x 2 / 3 . At 
T = 0, 2.5, 5, 7.5, the x m = 1, 1.0157, 1.0598 and 1.1251, while K = -42a; 2 / 3 + 0.705T 2 /a; 2 / 3 + 207a; has values 165, 
172, 193, 224 MeV, respectively. 

At fixed entropy, the second derivative of E(R)/A has a very different behavior than at fixed T. Namely, it 
decreases with R for a given S/A. This can easily be seen by noting that E X (R, S) = 4.836(S/ A) 2 A 2 ^ 3 / R 2 compared 
to E X (R,T) = 0.0517T 2 R 2 /A 2 / 3 . Here we used S = ^-TA/E F at low T. Because R now is in the denominator, the 
first derivative is negative and kq will decrease 



R 



2 d 2 (E(R)/A) 



dR 2 



= K (a) + 9.672(S/A) 2 A 2 / 3 /R 2 {Sk(a)/n (a) + 5} (15) 

s 

x m = (l + 9.672(S/A) 2 A 2 / 3 /i? 2 )" 3 (16) 



Note here that x m < 1, i.e., the saturation density is lower than p in an adiabatic process in contrast to an isothermal 
process. We have the following final results: 

re Q (a = 0) = 165 - 30(S7A) 2 
K Q (a = l/3) = 234 - 38(S/A) 2 

K Q {a = l) = 373 - 53(S/A) 2 (17) 

in MeV. 

At higher T, the nearly degenerate Fermi gas kinetic energy term is replaced by a virial expansion in pA 3 , where 
A is the quantum wavelength. Namely, Ek/A = (3/2)T(l + ^c„(pA 3 /4) n ) where the sum over n = 1, 2, 3, • • • has 
coefficients: ci = 1/2 5 / 2 = 0.177, c 2 = (1/8 - 2/3 5 / 2 ) = -3.3 x 10~ 3 , c 3 = 1.11 x 10" 4 , HQ. These coefficients 
come from antisymmetrization effects only. The role of clusters in the virial expansion is discussed in Ref.0. Since 
the c ra 's become small rapidly, we will keep terms up to C2. Then re is given by: 

k{T) = -(2196/T 2 )a4 + 9a(l + a)b a x 1 + a (18) 

The x m is again the minimum of E/A, but now evaluated with the new kinetic energy. The x m is affected by both 
the ci and ci terms. A limiting value of re can be obtained by taking T very large where an ideal gas Ek/A = (3/2)T 
leads to very simple results Q- 

We also investigate the case of constant entropy in the ideal gas limit using the Sackur- Tetrode law 
5/2 - ln(A 3 p/4). This law connects T to p or V as T = C s p 2 ' 3 . Here C s = [2n{hc) 2 /(mc 2 )} exp[(2/3)(S7A) - 5/3]. 
Then E K /A = (3/2)C S/ o 2/3 (l + p 2/3 E c n {[2-rth 2 /mC s ] 3/2 ) n )- The resulting E{R)/A has a structure similar to the 
result for a degenerate Fermi gas since both have a p 2 ! 3 dependence for the kinetic energy term but with different 
coefficients. This feature and a similar result at lower T suggests a duality in the energy per particle EOS at constant 
entropy and its associated hq and the T — EOS and its associated constant T re Q. 

It should be noted that the a; m 's for the constant T case are usually close to 1 except at high T, and greater then 
1, so that the minimum of the energy occurs at a density that is greater than po. This feature is not the case for 
the Helmholtz free energy F = E — TS. Since P = —dF/dV, F is minimum when P = 0. At low T, the entropy 
S — (it 2 /2)TA/ Ep, and thus Fk — Ek — TS, the part of the Helmholtz free energy coming from the kinetic energy 
part of E, is F K /A = (3E F /5)(l - (5tt 2 / 12) (T/E F ) 2 ) . The complete F/A = F K /A + U/A, where U/A is given by 
Eq.(J2Jl. The free energy per particle at low T can be obtained from Eq.@ by just changing the sign of the T 2 term: 
F/A = 21x 2/3 - (7r 2 /140)T 2 /a; 2/3 - b x + b a x 1+a . As with the E/A equation, this F/A equation is only valid when 
T <C Ef. For higher T, we use the virial expansion for Ek involving the coefficients c„ mentioned above. Similarly, 

the TS/A = (5/2)T - T \n(\ 3 p x / 4) + Cl {X 3 p x/A)/2 H H El. We note that the minimum of the free energy 

F(T, V) and energy E(T, V) at constant T don't coincide except at T = 0. However, the minimum of free energy 
F(T, V) at constant T and the minimum of energy E(S, V) at constant S are at P = 0. The minimum of F is at 
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subnuclear density, or x m < 1, at non zero T. For a = 1/3 and at T — 0, 2.5, 5, 7.5 MeV the x m = 1, 0.988, 
0.951, 0.878 and the curvature function « = 9x 2 n (d 2 (F/A)/dx 2 n ) = 234, 225, 199, 150 MeV, respectively. The F 
also has a maximum at even smaller x than its minimum point. The pressure P = pxd{F / A) / dx = p$x 2 d(F / A) / dx 
and is P = 2.1x 5 / 3 + 0.00705TV/ 3 - 0A5b o x 2 + 0.15(1 + a)b a x 2+a for T < £J F . At high T and/or small .t, the 
P = p a xT{l + J2c n (\ 3 Po/4) n x n }-0.15b o x 2 +0.15(1 + a)b a x 2+a . The spinodal line is obtained by setting dP/dV = 
or d 2 (F/A)/dV 2 = 0. The critical point occurs at the top of the spinodal line where d 2 P/dV 2 = also. Similar 
remarks also apply to E(S, V) with S held constant and for different values of S. The specific heat at constant volume 
is also easily obtained from the above results. At low T, Cy is linear in T and is (it 2 /2)T/ Ep. The behavior of Cy 
at higher T can be obtained from the virial expansion of Ek/A given above. At very high T the Cy — > 3/2. The 
T = limit where Cy = is connected to this very high T limit of 3/2 in a monotonic way with increasing T with 
no peak or cusp in Cy. This is to be compared with the strong peak found in Cy when inhomogeneities are included 
fl6|. In this latter case, cost functions arising from the surface energy associated with the formation of clusters are 
responsible for the very large peak as mentioned. These calculation [ly were done at a freeze out density of about 
1/3 nuclear matter density. At high density, i.e., above nuclear matter, heated matter may first form bubbles in a 
uniform background of liquid. The formation of a bubble has associated with it a surface energy. This aspect will be 
considered in a future paper. 



III. EFFECTIVE MASS EFFECTS 



An effective mass m* changes the kinetic energy term. If we parameterize it as m*/m = 1/(1 + r^p/po) 7 ), then 
at zero temperature the E/A will simply read E/A = (E K /A)x 2/3 + r(E K /A)x 2/3+ ' 1 - b x + b a x 1+a . The E K /A = 
(3/5)Ep = 21MeV, the value with the unrenormalized mass, and the effective mass correction appears in the additional 
second term with a 7 + 2/3 power dependence. The r is a constant and has a value of 1/2 when m* /m = 2/3 at po- 
If a = 7 — 1/3, then this additional effective mass term could be simply added to the b a term. Since the coefficients 
are determined to give the correct binding at p = p or x = 1 where E/A has a minimum, the b a coefficient must be 
adjusted also. The final incompressibility is still given by Eq.|0J) with = 7 — 1/3. If we parametrize m*/m with 
7 = 1 then for a = 2/3 the effective mass term can be added to the b a term and the final incompressibility is now 
given by Eq.Q with a — 2/3. Similarly if we parametrize m*/m with 7 = 2/3 then for a = 1/3 the effective mass 
term can be added to the b a term and the final incompressibility is now given by Eq.(@J) with a — 1/3. If 7 = 1/3, 
the r dependent effective mass term can be added to bo term and there is no effect on k and Sk. For general a and 
m*/m: ab a = (E B /A) + (1 - [3 7 - l]r)(E K /A)/3 and b = (1 + r){E K /A) + (E B /A) + b a and 

«o(a,r) = 9E B /A+(l + [3j-l} 2 r)EK/A + a(9EB/A + 3(l-[3j-l}r)E K /A) 

-» 3(55 + 69a + 7r [3 7 - 1] ([3 7 - 1] - 3a)) (19) 

For 7 = 1/3 or for a = 7 — 1/3 this k is r independent for reasons already given. When we parametrize m* /m with 
7 = 1, for a = 1/3, k = 234 + 42r and for r = 1/2, k = 255MeV. For a = 1, k = 372 - 42r -> 351MeV at r = 1/2. 
The skewness is also affected by the effective mass: 

Sk = -27(E B /A)(3 + 4a + a 2 ) - (E K /A)(U + 36a + 9a 2 + r[3 7 - 1]([3 7 - 1] - 3a)([3 7 + 11] + 3a)) 

-» -3(509 + 828a + 207a 2 ) - 21r[3 7 - 1]([3 7 - 1] - 3a)([3 7 + 11] + 3a) (20) 

For 7 = 1/3 or for a = 7 — 1/3 this is independent of r and for a < 7 — 1/3, Sk/ kq decreases with r; for a > 7 — 1/3, 
Sk/ 'kq increases with r. 

The effective mass also appears in the temperature dependent part of the energy per particle, modifying the 
additional term E X (R,T) to now be given by E X {R,T) = 0.0517(T 2 i? 2 /A 2 / 3 )(l/(l + r(R /R) 3 ' ( )). The last factor in 
this expression is the change the effective mass produces. The result of Eq.© can be used to calculate the new k. 
This leads to the following equation: 



k = Ka(a,r) 



-0.0705T 2 



(2 + ([3 7 - 3] [3 7 + 2] + 10)r- [3 7 - l][3 7 + 2]r 2 ) (2 + [3 7 + 2]r) f Sk(a,r) 



(1 + r) 3 (1+r) 2 \K (a,r 



(21) 



For the case of 7 = 1, at T = 7.5MeV and r = 1/2: k = 283MeV for a = 0, 335MeV for a = 1/3, 388MeV for 
a = 2/3 and 442MeV for a = 1. The T = 0, r = 1/2 values of n are 207, 255, 303, and 351 MeV and values 
of Sk are 2115, 2739, 3459, and 4275 MeV for each of these respective a values. For a fixed entropy, E X (R,T) — 



6 



A.836{S/A) 2 {A 2 / 3 /R 2 )(l/{l + r(R /R) 3 ^)) and thus 
hq = K (a,r) 

+3.548(S/A) 2 



(10- ([37 -2] [37 + 9] -2)r+ [3 7 - 2] [3 7 - 5] r 2 ) (2 - [3 7 - 2]r) fSk(a,r) 



(1+r) 3 (1 + r) 2 V«o(a,r) 



(22) 



At high T with effective mass, Er/A = (3/2)T(l + c„(p(l + r(p/p ) 7 ) 3/2 A 3 /4)™) due to the mass dependence 
in the De Broglie wavelength A. Due to the density dependence of the effective mass the c\ term contributes directry 
to incompressibility now. Then k is, up to c\ term, given by: 

<T) = rV2(i + T .^)i/2 ra ""2 7 + 7 + (1 + gT)^] + 9a(l + Q )6 a ^+ Q (23) 

The x m is again the minimum of E/A, but now evaluated with the new kinetic energy. At fixed entropy E^/A = 
(3/2)C s p 2/3 (l + X>„([27r?i 2 /mC,s] 3/2 (l + r(p/p ) 7 ) 3/2 )™)- The lowest order term having effective mass dependence 
is the c\ term. Up to this term, 

o /o »-2n 3/2 2/3 

3 „ 2/3 / 27Tft \ Cilm 



Kq( 5) = -3Csp 2 /3 xU 3 + \CspT ( 



mC 5 y (l+rx 7 ) 1 ^ 

~ 9 V 2 x^ 7 ) +9a(l + a)fo Q ^ +Q (24) 



(277 2 + 9 7 - 8) . ^ (97 + 4)(97- 2, 2 . ^ ^ 



IV. EXCLUDED VOLUME EFFECTS 

Excluded volume effects appear in the Van der Waals approach and modify the adiabatic incompressibility, leading 
to an equation of state of the form, for a given 5, 

E/A = (3/2)C s [p/(l - bp)} 2 / 3 - a aP (25) 

where C s = [2Tt{hc) 2 / {mc 2 )] exp[(2/3)(S/A) - 5/3] which connects T to p or V as T = C s/ 9 2/3 for ideal gas. An 
expansion of 1/(1 - bp) 2 / 3 = 1 + (2/3) (fy?) + (5/9)(&p) 2 + • • • leads to an infinite series of repulsive terms with various 
powers of the density. The numerator has p 2 / 3 so that additional terms such as p 5 / 3 , p s / 3 , p 11 / 3 , ■■■ are generated 
in this expansion. These terms collectively give rise to a much high incompressibility. If we had truncated the series 
at p 5 / 3 , then the resulting EOS would be similar to a Skyrme with a = 2/3. A Skyrme with a — 2/3 has an 
incompressibility of 303MeV. 

To explore higher density repulsive terms in a T = Skyrme parametrization, we can use the duality bewteen 
the EOS at constant entropy and the EOS at zero T 0. Thus we take the specific form of Eq. l|25|l with (3/2)Cg 
replaced by a constant parameter clk at zero T as a very simple model for studying excluded volume like effects in 
our E/A. We can fix the coefficient clq and an to give a binding energy per particle of 16MeV at p = 0.15/fm 3 : 

b K = a K (j-^) = 48 (sT^l) and b o = aoP = ( 3b f-i) ■ Thc resulting k is 

K = 96(66p-l)/((36p-l)(l-6p)) (26) 

This k has a minimum value of 754MeV at b = 3.60. Alternatively if we fix aK to be the T = Fermi kinetic 

2/3 

energy coefficient, cikPo = 21MeV, and determine ao and b by having correct binding energy at saturation we get 
b = 3.61 and an incompressibility of 754MeV. The two methods result in very similar values for b and k. A high 
incompressibility of 754MeV is equivalent to having a single Skyrme repulsive term with a = 2.85 from Eq.Q. This 
high incompressibility is unreasonable but it can be moderated by having higher order attractive terms in the density 
which, for example, can come from higher order attractive correlations of nucleons 0- The skewness is given by 



Sk = 192 



bp-1 
36p-l 



6 + 31(^0+45^-^1 +20^ bp 



1 — bp J \1 — bp J \ 1 — bp 



The effects of temperature and entropy dependences coming from a kinetic energy can be obtained from Eas. i|lU|l and 
(|15fl with the kq and Sk given above. 
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V. ASYMMETRIC AND FINITE SYSTEMS 

For N, Z asymmetric systems with Coulomb and finite size effects included, the T = K is reduced This was 

qualitatively explained in Ref.[3j as a reduction in the binding energy term in Eq.@. Detailed calculation done in 
Ref.[4| confirmed this feature. In this part of the paper we will present some very simple equations for the contribution 
of surface, symmetry and surface symmetry, and Coulomb and surface Coulomb. The surface, symmertry energy and 
Coulomb terms each reduce the binding energy per particle from 16MeV to 8MeV. For light nuclei around A = 27, 
the surface energy/^4 ~ 6MeV, the symmetry energy /A ~ OMeV and the Coulomb energy/A ~ 1.5MeV. For heavy 
nuclei around A = 216, the above numbers are changed to ~3MeV, ~lMeV and ~4MeV for the surface, symmetry 
and Coulomb terms, respectively. We will find a somewhat similar pattern of relative importance in k, but the surface 
term will be dominant even for large nuclei. The Coulomb is of next importance followed by the symmetry term for 
nuclei near the valley of beta stability. Surface symmetry and surface Coulomb terms will also be included. 



A. Surface Contribution 



First, let us investigate the role of finite size or surface correction, before considering the full complexity of the 
problem with Coulomb and symmetry energy terms included in the energy. The surface energy can be studied by 
using p(r) = p c /(exp[(r — R)/o] + 1) as a density profile. Here p c is the central density and a is a parameter that 
determines the surface thickness. This density can then be used in a Skyrme energy density functional: E(p) — 

&kP 5 ^ 3 — a-oP 2 + a a p 2+a - The clk is the kinetic energy constant and is taken as axp^ 3 = (3/b)Ep = 21MeV. The 
ao and a a are given above. The energy per particle is the integral of E(p) over all space divided by A. The resulting 
integrals can be simplified by an integration by parts and very accurately approximated by noting that the derivative 
of the Wood-Saxon profile gives a narrow peak of width a in the vicinity of R. A substitution y = (r — R)/a and 
a change in the lower limit of integration from —R/a to — oo gives a very accurate approximation to these integrals 
when R a, as is usually the case. As an example, the energy per particle from the ao attractive interaction 
term is (ao/^4)(47ri? 3 /3)p 2 [l — 3(a/R) + ir 2 (a/R) 2 — ir 2 (a/R) 3 }. The central density p c , R and a are also connected 
since the integral of p over all space must give A. Using the integration procedure just mentioned, the condition 
(47ri? 3 /3)p c (l + ir 2 (a/R) 2 ) — A is obtained. Note that in this last result the linear and cubic terms in a/R are absent 
in the constraint equation. This is not the case for the integrals involving the kinetic energy density, the ao and the 
a a terms. Each of these has a form involving surface thickness corrections that is (1 — ci(a/ R) + c 2 (a / R) 2 — Cy,(a/ R) 3 ). 
For the kinetic energy term c\ = 2.2769, C2 = 9.1046, C3 = 7.805. For the ao term c\ = 3, c 2 = C3 = tt 2 , and for the 
a a term, for a = 1/3, ci = 3.585, c 2 = 10.819 and c 3 = 11.642, for a = 1, c x = 9/2, c 2 = tt 2 + 3 and c 3 = 3n 2 /2. The 
energy per particle can be rewritten as: 

E(R)/A = 28.626C K (A 2 ^ 3 /R 2 ) - (58.89 + 36.61 / a)C (A/ R 3 ) + (23/a)C Q (A/(0.8565i?) 3 )" +1 (28) 

The Ck, Co and C a are the surface correction factors, each involving their corresponding ci, C2, C3 terms as [1 — 
a(a/R) + C2(a/R) 2 — C3(a/R) 3 ]/[1 + Tr 2 (a/R) 2 y , where the denominator arises from the constraint condition between 
R, p c and A. The j — 5/3 for the kinetic term, J = 2 for the ao term and j = a + 2 for the a a term. For a/R, <C 1, 
the surface energy can be approximated by keeping linear terms in a/R to give 

E diff{R)/A = -(a/R)(21 x 2.2769(i? /i?) 2 - (37 + 23 / a)3(R / R) 3 + (23/a)(3.585)(i? /i?) 3Q+3 ) 

= -(a/i?)(65.18A 2/3 /i? 2 - 506.11^/i? 3 + 459.65A 4/3 /i? 4 ) (29) 

for a = 1/3. This E diff -> as a -> and E diff = 10.36^ 2 / 3 MeV at a = 0.53fm and p = 0.15fm- 3 . To have 
correct surface energy of ~ 20A 2 / 3 MeV, we also consider a gradient dependent energy density of B(Vp) 2 = —BpV 2 p. 
With the Wood-Saxon density profile this becomes 

p c ll + (a/R) 2 (n 2 /3~2)} p 3 A 3B a A 

Ejfjrnd/ si — tj 77- : 777777 K3 £> = £> - 



2i?a [1 + 7r 2 (a/i?) 2 ] 2 2i?a 87ra R 4 %ira 2 R R 3 



3B 



/inpo\ 4/3 /JZo 



,2 



A 1 / 3 \ 3 J \R 



(30) 



The approximate form of this term is for a/R <C 1 and is the same form as the ao term in Ediff- The total surface 
energy then becomes 

E surf = 4nR 2 <j{p) = 4Trr 2 A 2 / 3 a{p) = E diff + E grad (31) 



8 



The surface energy has a contribution from the gradient term that varies with a as l/o while the diffuseness term 
varies directly with a. Note that as the surface gets sharper or a — > 0, the gradient term becomes very large and 
the decrease in binding from the surface diffuseness becomes small. Alternatively, as the surface thickness a becomes 
large, the gradient term is reduced while the other term increases. These behavior implies an extremum condition for 
E S urf versus a. We take this extremum to be at a = 0.53fm and use this to determine B. This extremum condition on 
E sur f leads to the result that the gradient term and the the reduction in binding energy from diffuseness make equal 
contributions to E sur f at its saturation. For a — 0.53fm, with the approximated form of E gra d, the B = 85.47MeV-fm 5 
and E di ff = E grad = 10.36A 2 / 3 MeV. The surface energy E sur f is the same form as Eg. (|2H with modified coefficient 
of ao term to be 542.43, i.e., 

E surf (R)/A = -(a/R)(65A8A 2/3 /R 2 - 542.43A/i? 3 + 459.65A 4/3 /i? 4 ) (32) 

The incompressibility will also have two contributions each from Ediff an d E grae i- To have the same surface energy 
Esurf = Ediff + Eg ra d for other values of a, a can be determined to have E d if / = — (47rpo/3) 1 ^ 3 (a/A 1 / 3 )(21 x 2.2769 — 
(37 + 23/a)3 + (22>/a)C a \) = 10.36/A 1 / 3 with C a \ the coefficient of a/R term in C a . Then the parameter B can be 
set by E grad = gi* (^) 4/3 ^ = 10.36/^/3. 

The values of n at the minimum E/A can be obtained from Eq.JJHI) with E gra d included. The energy E{R)/A with 

-Bgrad IS 

E(R)/A = 21C K (a){R /R) 2 /D 5 / 3 - (37 + 23/a)C {a)(R( ) / R) 3 / D 2 + {23/a)C a (a)(R / R) 3a+3 / D a+2 

+10MC G {a){R /R) 4 /D 2 (33) 

where 



C K {a) 


= 1 


- 2.2769(a/i?) + 9.1046(a/i?) 2 - 7.805(a/i?) 3 


Co (a) 


= 1 


- 3(a/R) + Tr 2 {a/R) 2 - 7r 2 {a/R) 3 


Ci/ 3 (a) 


= 1 


- 3.585(a/i?) + 10.819(a/i?) 2 - 11.642(a/i?) 3 


Ci(o) 


= 1 


- (9/2)(a/i?) + (tt 2 + 3)(a/i?) 2 - (37r 2 /2)(a/i?) 3 


C G (a) 


= 1 


+ (7r 2 /3-2)(a/i?) 2 


D(a) 


= 1 


+ n 2 (a/R) 2 



Since E sur f is not a small quantity we need exact calcuration here. For A = 216, 125 and 64 with a = 0.53fm, the 
respective values of k are 145, 128 and 103 in MeV, at the minimum values R m of E/A R m — 6.888, 5.709 and 4.527 
or x m = 1.0518, 1.0691 and 1.0981 respectively. These values are well below the value of k = 234 without finite size 
corrections. An approximate expression for these k is given by 

K = k(0) - 2.26k(0)/A 1 / 3 (35) 

while approximate calculation up to first order of a/R = (4:Trpo/3) 1 ^ 3 (Ro/R)(a/A 1 ^ 3 ) gives k — k(0)(1 — 2.01/A 1 / 3 ). 
The form of this result not only applies for a = 1/3, but also for softer EOS or smaller a, and harder EOS or larger 
a. The k(0) is n at R — > oo which is the infinite A k. Thus k(0) = 165, 234, 372 for a = 0, 1/3, 1 respectively. 
The simple result given by the equation comes from the observation that the surface correction coefficient scales with 
the k(0) coefficient for different a. This observation came from the detailed calculation that were done for various A 
from A = 64 to A = 216 and various values of a. For example, when a — 1 re(0) = 372 and detailed calculations gave 
k = 249, 224, 188 for A = 216, 125, 64 with a = 0.478fm which can be approximated as re = k(0)(1 - 1.98/A 1 / 3 ). 

The temperature and entropy dependence of the incompressibility can be evaluated using Eq. 10 for the T dependent 
kinetic energy with the kq and Sk determined above for the surface energy included. Beside the kinetic energy, the 
surface energy itself can have a temperature dependence which needs further investigation. 



B. Symmetry Energy Contribution 

We next investigate the role of symmetry energy terms. Symmetry energy terms appear in both the kinetic energy 
and potential energy parts of the total energy. The symmetry energy appears as (N — Z) 2 /A 2 in the binding energy 
per particle in the well known mass formulae. At low temperatures, T <c Ep, and for proton fractions y — Z J 'A, such 
that \2y — 1| <C 1, the kinetic energy per particle is 



E K /A = 2lx 2 / 3 {l + (5/9) (2y - l) 2 ) + (7r 2 /140)(T 2 /a; 2 / 3 )(l - (l/9)(2y - l) 2 ). 



(36) 
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At higher temperatures, an expansion in pA 3 is more appropriate and the kinetic energy per particle now becomes 


E K /A = (3T/2){1 + c 1 (A 3 jO0 /4) a ;(l + (2y - l) 2 ) + c 2 (A 3 p /4) 2 x 2 (l + {2y - if) + ■••}. (37) 

The y dependence in Ea. (|37J) . to the expansion order c 2 , is not restricted to \2y — 1| <C 1, i.e, it is valid for any y. 
The potential energy per particle, U/A, is [l4T | 

U/A = -b n x + b a x 1+a + (2/3)(l/2 + x )b Q (2y - ifx - (2/3)(l/2 + x 3 )b a (2y - l) 2 x a+1 (38) 

The xq is a parameter which we fix to give the empirical symmetry energy coefficient in the binding energy mass 
formula. This symmetry energy coefficient is taken to be 25 MeV. Note that the kinetic energy contributes to this 
value and from Ea. (|36J) is (35/3) MeV at x = 1 and T = 0. The remaining (40/3) MeV comes from the third and 
fourth term on the rhs of Eq.lj^J). Thus (40/3) = (2/3)(f/2 + x )6 - (2/3)(l/2 + x 3 )b a then x and x 3 satisfy 
(37 + 23/o)xq — (23/a)x3 = 3/2. For a = 1/3, 106xo — 69^3 = 3/2. If x 3 — I, which is the value taken for 
Skyrme forces ZRI, ZR2, ZR3, then xq = 0.6651. If x 3 = —1/2, so that the fourth term in Ea. l|38|) is absent, then 
Xq = — 0.3II3. If Xq = —1/2, so that the third term is now absent, then x 3 = —0.7899. The nuclear incompressibility 
involves the second derivative of E/A with respect to x or p so that the fourth term in Ea. (|38|l . which involves x a+1 , 
is very important in the asymmetric part of the nuclear incompressibility. Terms involving linear powers in x do not 
make explicit contributions, but they do affect the minimum point in the evaluation of E/A. However because of the 
uncertainty in the value of x 3 , and thus in the coefficient 1/2 + x 3 that appears in the fourth term, the contribution of 
the asymmetry term to the nuclear incompressibility would seem to be somewhat uncertain. This turns out not to be 
the case because of cancellations so that the coefficient multiplying the (1/2 + x 3 ) part in k is small as will be found. 
We can evaluate the relative importance of the contribution of the symmetry energy to k by isolating those terms 
associated with it and turning off the Coulomb and surface contributions. Let us consider the choice of a general x 3 , 
which has an associated xq to give correct symmetry energy. Our E(R) with symmetry terms included and is 

E{R)/A = C K (y)(R /R) 2 - C (y)(R /R) 3 + C a (y)(R /R) 3{a+1) (39) 

with 

C K (y) = 21 + (35/3)(2y-l) 2 

C (y) = (37 + 23/a)- (40/3 + (46/3a)(l/2 + x 3 ))(2y - I) 2 

C a (y) = (23/ a )-(46/3a)(l/2 + x 3 )(2y-l) 2 (40) 

Thus 

= -2C K (y)(R /R) 2 + 3C (y)(Ro/R) 3 - (3a + 3)C a (y)(R /R) 3a+3 (41) 
K = 6C K (y)(R /R) 2 -12C (y)(R Q /R) 3 + (3a + 3)(3a + A)C a (y)(Ro/R) 3a+3 (42) 
Sk = ~24C K (y)(R /R) 2 + 60C Q (y)(R Q /Rf - (3a + 3)(3a + 4)(3a + 5)C a (y)(R /R) 3a+3 (43) 

For a ~ 1/3, the minimum point of E{R)/A is determined by a simple quadratic equation. The value of n at this 
minimum is accurately given by, up to (2y — l) 2 order, 

k = 234 - (300 - 16(1/2 + x 3 ))(2y - l) 2 (44) 

withx TO = (R /R m ) 3 = (l-(0.271-0.197(l/2 + a;3))(2 2 /-l) 2 ) 3 . The small coefficient of 16 compared to 300 in front 
of the term involving x 3 , means that uncertainties in x 3 do not have a pronounced effect on the value of n. Similar 
results are also true for other values of a, with the second 1/2 + x 3 term about 5% of the first term in the (2y — l) 2 
dependent part of k. For general a, using either Eci. (|4*2|) or ©, 



_ 182.9 + 403.33. + 190a 2 _ (1/2 + x 3 )(9.333 + 28a) _ 

0.7971 + a v y ' 0.7971 + a 



with x m = (1 - (0.3060- 0.2222(1/2 + x 3 ))(2y- l) 2 /(0. 7971 + a)) 3 . Surface symmetry corrections can also be easily 
included following the procedure given above. Since the volume and volume symmetry terms have the same density 
dependence, the surface corrections to each will scale in exactly the same way with their corresponding volume term. 
Thus we can rewrite the symmetry term in the above equation as (300 — 16(1/2 + x 3 )){2y — 1) 2 (1 — 2.26) / 'A 1 / 3 ) for 
a = 1/3, where the 2.26/A 1 / 3 part of this result is now the surface symmetry correction to the symmetry energy (c.f, 
Eq.^SJ). 
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At non zero but low T, Eq.J2U) gives a small T dependent kinetic energy term E x = (it 2 /140)T 2 (i?/i? ) 2 (l — (2y 
l) 2 /9). Using Eq.©, 



k(T) = K (0) - (tt 2 /70)T 2 (1 - (2i/ - l) 2 /9)(S , fc/ K (0) + 1) 



(46) 



where k(0) is the incompressibility at zero T for an asymmetric system which is given by Ea. (|45|l . On the other hand 



E x = (35/tt 2 )(S/A) 2 (R /R) 2 (1 
incompressibility of 



(2y — 1) /9) for a given entropy since S/A = (n /2)T/Ep . This gives the adiabatic 



kq(T) = k(0) + (70/7r 2 )(5/A) 2 (l - (2y - l) 2 /9)(Sfc/«(0) + 5) 



(47) 



At higher T, the energy is the sum of Eq. (|37|l and Eq. (|38|l and we can neglect higher order terms in virial expansion 
keeping terms only up to c-i. Then k is given by: 



k(T) = 27 



C2 ( 27T?l Z 
2^2 



(f ) 2 xl l (l + (2y-l) 2 ) + 9a(l + a)b a 



1 



+ x 3 (2y-iy 



(48) 



The x m is again the minimum of E/A, but now evaluated with the new kinetic energy. The x m is affected by both 
the ci and ci terms. A limiting value of k can be obtained by taking T very large where an ideal gas Ek/A = (3/2)T 
leads to very simple results. Adding this Ek/A to U/A given by Eq. l|38|) leads to a minimum x m given by x^ = 
b [l - (2/3)(l/2 + x )(2y- 1) 2 ]/(6 Q (1 + a)[l - (2/3) (1/2+ #3) (2y- l) 2 ]). In this high T limit k is given by the second 
term on the right side of Eq. I|48|l and is 



«(T ~ lOMeV) = 9ab 



i-¥ 

3 V 2 



x )(2y-lf 




(49) 



At fixed entropy E K /A = (3/2)C s p 2/3 (l + J2 c n([27rh 2 /mC s } 3/2 (l + (2y - l) 2 )). The resulting E{R)/A has a struc- 
ture similar to the result for a degenerate Fermi gas since both have a p 2 ^ 3 dependence for the kinetic energy term 
but with different coefficients. This feature and a similar result at lower T suggests a duality in the energy per 
particle EOS at constant entropy and its associated kq and the T = EOS and its associated constant T K even for 
asymmetric system. 



C. Coulomb Energy Contribution 

Next, let us consider the Coulomb energy contribution. For a uniformly charged sphere we have E c /A = 
(3/5)e 2 (Z 2 /A)/R = (3/5)(e 2 /hc)(hc)(Z 2 /A)/R. As before, we turn off the other two contributions, now the sur- 
face and symmetry, to first study the isolated importance of the Coulomb interaction. In this case our energy per 
particle at T = is 

E(R)/A = 28.627A 2/3 /R 2 - (58.89 + 36M/a)A/R 3 + (23/a)(A/(0.8565 J R) 3 ) Q+1 + 0M28(Z 2 /A)/R (50) 

The Coulomb contribution to k can be expressed in terms of a simple result, using Eq.JfjJ, 

k c = 4E C /A + (E c /A)(Sk/ Ko ) (51) 

The skweness Sk, third derivative of Eq(R), is negative. Here Eq(R) is the energy without the Coulomb correction 
and Ko its associated incompressibility. The rhs is evaluated at the minimum without Coulomb corrections. For a 
uniform sphere of radius R, for a = 1/3, the 

Hi = 234 — 4.75Z 2 /A 4/3 (52) 

For A = 216, Z = 86, R m = 7.14 and k = 207MeV; for A = 125, Z = 50, R m = 5.914 and k = 215MeV; for A = 64, 
Z=28, R m = 4.717 and k = 220MeV. If the density profile is taken to that of a Wood-Saxon form, the Coulomb 
energy can be approximated by (3/5)(Z 2 //l)(l — (7tt 2 /6)(a/R) 2 )/ R [l7) |. Thus surface contributions appear both 
directly and indirectly, modifying the other terms. Adding surface contributions gives, using Eq.JSJ for a = 1/3, 

k = 234 - 4.75Z 2 /A 4/3 + 22.9Z 2 /A 2 (53) 

Other values of a can also be easily evaluated. Neglecting surface-Coulomb we have for a = 1, k = 372 - 6.25Z 2 /A 4 / 3 
and for a = 0, K = 165 — 3.875Z 2 /v4 4 / 3 . The ratio of the Coulomb coefficient of Z 2 /A 4 / 3 to volume coefficient is not 
the same for each a so that a simple scaling result is not possible. 

Since both the Coulomb energy and the T dependent part of kinetic enrgy arc small the effects of these energies to 
the incompressibility are additive. 
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VI. CONCLUSIONS 



In this paper we investigated the behavior of the asymmetric nuclear matter incompressibility at finite temperature 
and entropy uing a mean field theory. The incompressibility is related to the curvature of the energy as a function of 
density, with the isothermal incompressibility having temperature held constant and the adiabatic incompressibility 
has entropy held constant. These two incompressibilities are the same at T = or S = but have very different 
behaviors with temperature T and entropy S. At all T or S these incompressibilities are very sensitive to the choice 
of a, the density dependent power of the repulsive interaction. We look at extreme cases of a very soft equation 
of state where a goes to zero to a very hard equation of state with a = 1. The range of k is from 165MeV to 
372MeV. For a moderate equation of state that is frequently used, a = 1/3 and k = 234MeV. When a temperature 
dependence is included, the value of the isothermal incompressibility k increases sharply with T. At low T or T 
much less than the Fermi energy, n increases as T 2 . For example at T = 7.5MeV and a = 1/3, k is 302MeV. By 
contrast, the adiabatic incompressibility was shown to decrease with increasing entropy and it eventually goes to zero. 
The adiabatic incompressibility is shown to arise from an equation of state (EOS) or energy per particle that has a 
structure that is similar to a T — Fermi gas. 

The incompressibility is also studied in finite systems where we investigated the role of surface, symmetry energy 
and Coulomb energy effects. First we studied each of these terms separately to illustrate the relative importance 
of them. Of these three terms, the surface had the largest effect, followed by the Coulomb and then the symmetry 
term for nuclei in the valley of (3 stability. Uncertainties in the X3 parameter of the Skrymc interaction were shown 
not to be so important. In our study of the surface term we used a Wood-Saxon density distribution. We found an 
approximation scheme to do various integrals (similar to the Sommerfeld method for Fermi gas integrals) that appear 
in our analysis and obtained very accurate analytic results for the surface contribution to the energy per nucleon. 

This work was supported in part by the DOE Grant No. DE-FG02-96ER-40987 and DE-FG02-95ER-40940 and 
in part by Grant No. R05-2001-000-00097-0 from the Basic Research Program of the Korea Science and Engineering 
Foundation. 
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